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The computation of radiative energy loss in a dynamically screened QCD medium is a key ingre- 
dient for obtaining reliable predictions for jet quenching in ultra-relativistic heavy ion collisions. We 
calculate, to first order in the opacity, the energy loss suffered by a heavy quark traveling through 
an infinite and time-independent QCD medium and show that the result for a dynamical medium 
is almost twice that obtained previously for a medium consisting of randomly distributed static 
scattering centers. A quantitative description of jet suppression in RHIC and LHC experiments 
thus must correctly account for the dynamics of the medium's constituents. 
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I. INTRODUCTION 



, Studying the suppression pattern of high transverse momentum hadrons is a powerful tool to map out the density 
' of a QCD plasma created in ultra-relativistic heavy ion collisions [1-3]. This suppression (called jet quenching) results 
from the energy loss of high energy partons moving through the plasma [4-7]. Recent non-photonic single electron 
data [8, 9] (which present an indirect probe of heavy quark energy loss) showed that radiative energy loss alone can 
not explain the results as long as realistic parameter values are used [10]. Inclusion of collisional energy loss [11-14] 
. improves the agreement with available data [15], but still does not yield a perfect description. 

The currently available studies suffer from one crucial drawback: The medium induced radiative energy loss is 
computed in a QCD medium consisting of randomly distributed but static scattering centers ( "static QCD medium" ) . 
In such a medium the collisional energy loss is exactly zero. This approximation was motivated by early estimates [16- 
21], which indicated that the typical collisional energy loss should be small compared to the radiative one. However, 
recent calculations [11-14] showed that the collisional contribution is important and comparable to the radiative 
energy loss. The static approximation is thus qualitatively wrong as far as the computation of collisional energy loss 
is concerned and should therefore also be revisited in the context of radiative energy loss. 

In this paper, we report on a first important step, the calculation of heavy quark radiative energy loss in an infinite 
. and time-independent QCD medium consisting of dynamical constituents. By comparing with the static medium 
calculation, this permits us to qualitatively assess the importance of dynamical effects on radiative energy loss. The 
more demanding problem of including finite medium size corrections, i.e. the Landau-Pomeranchuk-Migdal (LPM) 
effect, will be left to a future study. 

Here is the outline of our paper: In Section II we compute, to first order in the opacity, the radiative energy loss in 
an infinite, dynamical QCD medium. In Section III we obtain the corresponding result in the static approximation. 
While the analytical results in both cases lead to formally very similar expressions, they give remarkably different 
numerical values for the energy loss. These are presented in Section IV. We will see that a dynamical medium leads 
to approximately twice the radiative energy loss obtained in the static approximation. In Section V we present a short 
summary and conclude that representing the QCD medium by a random ensemble of static scattering centers is not 
a good approximation for RHIC and LHC phenomenology Some technical steps of our calculation are reproduced 
in the Appendix. We will use the following notation for 4-vectors: k = (fcrj,k) = (ko, k z , k), i.e. k (with an explicit 
vector superscript) describes a 3-vector while k (without a vector superscript) denotes the 2-vector transverse to the 
direction of motion z of the heavy quark. Correspondingly d 3 k = dk z d 2 k. 



II. RADIATIVE ENERGY LOSS IN A DYNAMICAL QCD MEDIUM 

In this Section we compute the medium induced radiative energy loss for a heavy quark to first order in the 
opacity. For simplicity we consider an infinite QCD medium and assume that the on-shell heavy quark is produced 
at time xq = — oo. In this medium we compute the radiative energy loss per unit length, For phenomenological 
applications in heavy-ion collisions one would, as a first approximation, use this result and simply multiply it with 
the effective thickness L of the medium to calculate the total energy loss. A more rigorous derivation would have to 
start with a finite size medium from the beginning; we leave this for the future. 

Medium induced radiative energy loss is caused by the radiation of one or more gluons induced by collisional 
interactions between the quark of interest and partons in the medium. The energy loss rate can be expanded in 
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the number of scattering events suffered by the heavy quark which is equivalent to an expansion in powers of the 
opacity. For a finite medium, the opacity is given by the product of the density of the medium with the scattering 
cross section, integrated along the path of the heavy quark. The lowest (first) order contribution corresponds to one 
collisional interaction with the medium, accompanied by emission of a single gluon. We adopt this as a definition of 
the "first order in opacity" also for the infinite medium. 

For a medium consisting of dynamical quarks and gluons in thermal equilibrium, the corresponding energy loss 
contribution involves two cut Hard-Thermal Loop (HTL) gluon propagators. The associated Feynman diagrams are 
plotted in Figs. 10-12 and computed in Appendices B-D. The diagrams represent an on-shell heavy quark with 
momentum p' which (in arbitrary order) exchanges a virtual gluon of momentum q with a parton in the medium and 
radiates a gluon with momentum k. The heavy quark emerges with (measured) momentum p. Since the exchanged 
gluon momentum is space-like [13, 18, 22]), only the Landau damping contribution (q < |q|) to the cut HTL effective 
gluon propagator D(q) needs to be taken into account [13, 17, 18]. 

The radiated gluon has timelike momentum k = (to, k), so only the quasi-particle contribution at uj > |k| from the 
cut gluon propagator D{k) contributes [23-25]. Energy and momentum conservation requires p' = p + k + q. Since 
our focus is on heavy quarks with mass M 3> gT, we neglect the thermal shift of the heavy quark mass. 

The effective gluon propagator has both transverse and longitudinal contributions [26-32]. The 1-HTL gluon 
propagator has the form 

where I = (?o, 1) is the 4- momentum of the gluon and P^vil) and Q^vQ.) are the transverse and longitudinal projectors, 
respectively. The transverse and longitudinal HTL gluon self energies Ht and n^, are given by [29] 
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where y = and /i = gTy/N c /3 + Nf/6 is the Debye screening mass. 

While the results obtained in this paper are gauge invariant [22], we present the calculation for simplicity in Coulomb 
gauge. In this gauge the only nonzero terms in the transverse and longitudinal projectors are 

P ij (i) = s ij + ^, Q oo (0 = -I = i - | = i - y 2 - (2.3) 

As in [33-39], we assume validity of the soft gluon and soft rescattering approximations (see Appendix A for details). 
With these assumptions we compute in Appendices B-D the diagrams M 10 , M^i and M lj2 , which contribute to the 
first order radiative energy loss. For the interaction rate we find 

r(E) = ^2ImM tot = ^(2ImM 1 ,o + 2ImM 1 , 1 + 2ImM 1 , 2 ), (2.4) 



where (see Eqs. (B21), (C10), and (D6)) 

2ImM 1;0 = 8Eg 4 T [t a ,t c ] [t c ,t a ] 
2ImM M = 8Eg 4 T [t a ,t c ] [t c ,t a ] 
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Here [t a ,t c ] is a color commutator, m 2 = yU 2 /2 is the effective mass for gluons with hard momenta k > T, and 
X = M 2 x 2 + m 2 where x is the longitudinal momentum fraction of the heavy quark carried away by the emitted 
gluon. We assume constant coupling g. 

By using the above equations, the interaction rate becomes 
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where we used [t a ,t c ] [t c ,t a ] = C2{G)CrDr (with 02(G) = 3, Cr = §, and Dr = 3) and, in the second step, the soft 
rescattering approximation |fc| <C k z w u>. 

The interaction rate sums over all initial and final colors of the heavy quark. The heavy quark radiative energy 
loss per unit length is obtained from the above expression for the interaction rate by weighting it with the energy lo 
of the emitted gluon and averaging over the initial color of the heavy quark [18, 22, 25]: 
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where the second step is obtained after angular integration and we defined a "dynamical mean free path" Ad yn through 

X^ n = C 2 (G)a s T =3a s T. (2.10) 

Under the assumption that a s is not running, Eq. (2.9) can be further analytically integrated over ^ I k | ^ fc max 
where fc max — 2Ey / x(l — x) [37]. We obtain 
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Alternatively, Eq. (2.9) can be integrated over < |q| < g max where <7 max = V '4ET [14]. This leads to an analytical 
expression for radiated gluon spectrum J7d yn (fc, %)■ We find 
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with 

G( X , k, k 2 ) = (fc 2 + X ) 2 - 2« 2 (fc 2 - X ) + « 4 (2.15) 

It is worth noting that each of the three contributions in Eq. (2.5) diverges logarithmically in the limit of zero transverse 
momentum of the exchanged gluon, q^O. The reason is that in a dynamical QCD medium both transverse and 
longitudinal gluon exchange contribute to the radiative energy loss. While Debye screening renders the longitudinal 
gluon exchange infrared finite, transverse gluon exchange causes a well-known logarithmic singularity [25, 40], due to 
the absence of a magnetic screening [41]. However, the infrared divergences cancel in the sum (2.4), giving rise to a 
finite energy loss rate. (This can be seen from Eq. (2.14), where analytical integration over q is performed.) This is 
a nontrivial result since it was previously believed that in a dynamical medium a magnetic cutoff must be artificially 
introduced in order to avoid divergent results [25, 40]. 



III. RADIATIVE ENERGY LOSS IN A STATIC QCD MEDIUM 



Let us now briefly revisit for comparison the heavy quark radiative energy loss in a static QCD medium, using 
a derivation which parallels that of the previous section and thus clearly exhibits the differences between the two 
situations. We again consider an on-shell heavy quark produced in the remote past propagating through an infinite 
QCD medium that now consists of randomly distributed static scattering centers. The static interactions are modeled 
as Gyulassy-Wang [33] color-screened Yukawa potentials 

V n = V(q n )e l ^ x " (3.1) 
= 27^X40 e-^"-" T an (R)®T an {n) , 

here x n is the location of the n th scattering center, the two T symbols (with a n being summed over) denote the color 
matrices of the heavy quark and the n th scattering center, and w(q„) = 47ra s /(q 2 + yu 2 ). 

The diagrams contributing to the radiative energy loss at first order in opacity are shown in Fig. 1. As seen in 
the figure, the quark scatters of one of the color centers with momentum q = (0,q z ,q) and radiates a gluon with 
momentum k = (w, k z ,k). As in the previous section, energy-momentum conservation requires p' = p + k + q. 



M stat 




Jq- 




§ 

FIG. 1: Feynman diagrams Mf tat , M| tat and M| tat contributing to the soft gluon radiation amplitude in a static medium to 
first order in opacity. The static color charge has color a and exchanges momentum q = (0,q z ,q) with the heavy quark. The 
radiated gluon has color c and carries momentum k = (uj,k z ,k). 

The procedure for the calculation of the diagrams shown in Fig. 1 was already presented in [37] (see in particular 
appendices A and B there), so we will not repeat it here. We obtain 
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where e = e(k) — ^0, x Q, e ^ iAj - j , ej (with e = (1,0) or (0,1)) is the polarization vector of the emitted gluon. By using 
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Eqs. (3.1), (A14), and (C7), together with p ■ e w p' ■ e w e • fe/x [34, 37], Eq. (3.2) becomes 
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Squaring this and ensemble-averaging the result over the positions Xi of the scattering centers gives 
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In the last step we used [34] Tr (T ai T a2 ) — S aia2 C2(T)d,T / d g (we assume that all scattering centers ("target partons") 
have the same dy-dimensional color representation with Casimir C^iT)). Furthermore, we assumed (similarly to [34]) 
that the ensemble average over the phase factor gives 
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where V = LA±, with L and A± being the depth (along the heavy quark's path) and transverse area of the medium. 

We can now compute the interaction rate for the heavy quark in a static medium of scatterers with color represen- 
tation T: 



rstat 
T 



(*)=/ 



= /(^^(^ 



k 2 +M 2 x 2 +mt 



2xE 



)^«f| 2 ) 



Aa s [t c ,t a ][t a ,t c }C 2 (T) 



dr 1 f dx d 2 k d 2 q 



d g V J x (2tt) 2 (2tt) 2 



b(0,q)| 2 



k+q 



X (k + q) 2 + X 



_ , Cj?a s 1 f dx d 2 k d 2 q 
D R d T — / 

TT At J X TT 7T 



fc 



fc+<? 



(g 2 + /i 2 ) 2 \k 2 + X (k+q) 2 + X 



(3.7) 



where we defined the mean free path for a heavy quark scattering off quark- type ( a q" ) or gluon-type ( "</" ) scattering 
centers through [15, 33, 34] 
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After averaging over the initial colors of the heavy quark and the scattering centers and weighting the rate with the 
energy u> of the emitted gluon, the heavy quark energy loss in an infinite static QCD medium is then given by 
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This leads to 
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where c(rif) = 6 1 ^" 2 pp^jyg is a slowly increasing function of n/ that varies between c(0) w 0.73 and c(oo) w 1.09. 
For a typical value n/ = 2.5 (which we use in this paper), c(2.5) s=s 0.84 ~ 1. 

As in the previous section, under the assumption that as is not running, Eq. (3.10) can be analytically integrated 
over k or q. Integration over k yields 
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with Q given by Eq. (2.15). Note that the emitted gluon spectra, Eq. (2.14) for Jd yn (k, x) and Eq. (3.15) for J s tat{k, x), 
while showing some similarities, don't exhibit a similarly simple analytical relationship as was the case for J7dyn(<7, x) 
and J sta ,t(q,x) (see Eq. (3.13)). 

Finally, we can compare the radiative energy loss rates to first order in opacity in dynamic and static QCD media, 
Eqs. (2.8) and (3.10). Both expressions yield an energy loss that increases linearly with the path length L through the 
medium. This reflects our neglect of LPM interference effects [4, 5] - our result corresponds to the Bethe-Heitler limit. 
In spite of this and many other similarities between Eqs. (2.8) and (3.10), they feature two important differences: The 
first is an (9(15%) decrease in the mean free path 
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which increases the energy loss rate in the dynamical medium by O(20%). The second difference is a change in the 
shape and normalization of the emitted gluon spectrum, which in the energy loss rate is reflected by the replacement 
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As we will see in the next section, this second difference leads to an additional significant increase of the heavy quark 
energy loss rate and of the emitted gluon radiation spectrum by about 50% for the dynamical QCD medium. 
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IV. NUMERICAL RESULTS 



In this section we present numerical results for the total radiative energy loss, as well as its differential rate with 
respect to the energy fraction and transverse momentum carried by the exchanged and emitted gluons, to first 
order in opacity, using the expressions for infinite QCD media derived in the two preceding sections. Specifically, 
we discuss a static quark-gluon plasma of temperature T = 225MeV, with nf = 2.5 effective light quark flavors and 
strong interaction strength a s = 0.3, as representative of average conditions encountered in Au+Au collisions at RHIC. 
Further below we will raise the temperature of the medium to T = 400 MeV to simulate (average) conditions in Pb+Pb 
collisions at the LHC. For the light quarks we assume that their mass is dominated by the thermal mass 
where fi = gT^/l+Nf/6 w 0.5 GcV is the Debye screening mass. The charm mass is taken to be M = 1.2 GeV, and for 
the bottom mass we use M = 4.75 GeV. As noted before, the radiative energy loss in the Bethe-Heitler limit considered 
here increases linearly with the path length L traveled by the fast quark - for normalization purposes we will set this 
path length to a standard value of L = 5fm throughout. 
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FIG. 2: Fractional radiative energy loss for an assumed path length L = 5 mi as a function of momentum for light, charm 
and bottom quarks (left, center, and right panels, respectively). Full and dashed curves correspond to the energy loss in a 
dynamical and a static QCD medium, respectively, and are obtained from Eqs. (2.8) and Eq. (3.10). 



In Figure 2 we compare the momentum dependence of the radiative energy loss over an assumed path length 
L = 5 fm for a dynamical and a static QCD medium (as given by Eqs. (2.8) and (3.10)) . For all three types of quarks, 
the dynamical medium is seen to cause almost 70% higher energy loss than the static medium. The left panel in Fig. 4 
below shows that <~ 50% of this increase arises from the larger values of the function Jdyn(k,x) which describes the 
shape of the emitted gluon spectrum, with an additional ~ 20% increase stemming from the shorter mean free path 
in the dynamical medium. 

To better understand the kinematic distribution of the dynamical medium effects we will now study the energy loss 
differentially as a function of the gluon energy fraction x = uo/E and the transverse momenta k and q of the emitted 
and exchanged gluons. We define the double-differential energy loss spectra 
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In Figure 3 we plot the transverse momentum dependence of the exchanged (S(q), left panel) and emitted gluon 
spectrum (S(k), right panel). The largest differences between static (dashed) and dynamical media (solid) are observed 
at low transverse momenta |q|, |fc| < 1 GeV, especially for the exchanged momentum spectrum S(q). The latter is 
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qJGeV] MGeV] 

FIG. 3: Transverse momentum spectrum of the exchanged (Eq. (4.3), left) and emitted (Eq. (4.4), right) gluons for charm 
quarks traveling for L = 5 fm through a medium of dynamical (solid) or static (dashed) scatterers. The initial charm quark 
energy was assumed to be 20 GeV. Note that \q\ and |fc| are denoted as qx and kx in the figure. 

seen to exhibit qualitatively different low-q behavior in static and dynamical media, but the difference is seen to 
rapidly disappear at |q| > 1 GeV (middle panel in Fig. 4), as expected from Eq. (3.17). 

While the dynamical medium effects on the exchanged gluon spectrum are mostly concentrated at low transverse 
momenta, the right panels in Fig. 3 and (more clearly) Fig. 4 show that the effect on the emitted gluon spectrum 
extends over the entire transverse momentum region, causing enhancements by more than a factor 2.5 at low |fc|, and 
settling down at high |fc| > 1 GeV to an approximately |fc [-independent enhancement by a factor <~ 1.3. 




FIG. 4: Left panel: Ratio of dynamical and static energy losses as a function of initial quark momentum p. The dashed 
line shows the 20% increase in the energy loss due to the shorter mean free path Ad yn < A s tat in the dynamical medium (see 
Eq. (3.16)). The dashed region indicates the additional ~ 50% increase in the dynamical energy loss due to Eq. (3.17). Middle 
panel: \q\ dependence of the ratio of the energy loss integrands J{q) = J dxj(q,x) in Eqs. (2.11), (3.12). Right panel: \k\ 
dependence of the ratio of the energy loss integrands J(k) = J dxj(k,x) in Eqs. (2.13), (3.14). Note that \q\ and |fe| are 
denoted in the figure as qx and kx respectively. 

The fc-integrated effect on the radiative charm energy loss is shown in the left panel of Fig. 4. The energy loss ratio 
between dynamical and static media is almost independent of the momentum p of the fast charm quark, saturating at 
~ 1.7 above p > 20 GeV and being even somewhat larger at smaller momenta. (This includes the ~ 20% effect arising 
from the shorter mean free path in the dynamical medium.) We checked that the dynamical enhancement persists at 
constant level to the largest possible charm quark energies. This is shown explicitly in the left panel of Fig. 5 where we 
plot the charm quark energy loss ratio for an otherwise identical medium of higher temperature T = 400 MeV ( "LHC 
conditions" ) . There is no quark energy domain where the assumption of static scatterers in the medium becomes a 
valid approximation. 

The mass of the fast quark plays only a minor role for its energy loss. The right panel in Fig. 5 shows (for LHC 
conditions, but similar statements apply at lower medium temperature) the asymptotic energy loss ratio for very high 
energy quarks as a function of the quark mass. While the dynamical enhancement is largest for light quarks, the 
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FIG. 5: Radiative energy loss in a medium of temperature T = 400 MeV ("LHC conditions"). Left panel: Ratio of the fractional 
radiative energy loss in dynamical and static media for charm quarks as a function of initial quark energy E. The ratio saturates 
quickly to a constant value above E ~ 100 GeV. Right panel: Asymptotic value of the radiative energy loss ratio for high energy 
quarks as a function of their mass, with marks indicating the light, charm and bottom quarks. 



difference between light and bottom quarks is only about 15%, and b quarks still suffer about 70% more energy loss 
in a dynamical medium than in one with static scattering centers. 

Stronger quark mass effects are seen in the shapes of the exchanged and emitted gluon spectra themselves. Figure 6 
shows the fractional energy loss as a function of the energy fraction x = lo / E of the radiated gluon relative to the 
initial quark energy. For both static and dynamic media the fractional energy loss for bottom quarks (i.e. the emitted 
gluon energy spectrum) is seen to be significantly softer for bottom than for charm quarks. No such strong mass effect 




FIG. 6: Differential fractional radiative energy loss S(x) as a function of radiated energy fraction x = uj/E to first order in 
opacity, for charm (left) and bottom quarks (right) of initial energy E = 20 GeV traveling for L = 5fm through a dynamical 
(solid lines) or static medium (dashed lines) of temperature T = 225 MeV. 



is visible in the shapes of the transverse momentum spectra of exchanged and emitted gluons: Figs. 7 and 8 show 
that, at fixed x, the exchanged and emitted gluon transverse momentum spectra have very similar shapes for charm 
and bottom quarks, and that the main difference shows up in the x-dependence of these spectra. This strong quark 
mass effect on the emitted gluon energy spectrum is a consequence of the well known "dead-cone effect" [37, 43]. 
Except for very low x- values, the differential energy loss S{x) is a decreasing function of \ = M 2 x 2 +m 2 ,, and \ grows 
significantly faster with x for bottom than for charm quarks. Thus, as x increases, the contribution to the energy loss 
decreases more rapidly for bottom than for charm quarks. This is borne out by Fig. 6. 

In the present study, large- a; contributions to the total energy loss are not strongly suppressed, and even for bottom 
quarks the contribution from x regions where the soft gluon approximation ui <C E becomes doubtful could be as large 
as 30%. This is mainly a deficiency of our approximations - for a medium of static scatterers it is known that the 
LPM effect strongly suppresses the emission of large- x gluons (see Eq. (11) in the second paper of Ref. [37]). Including 
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FIG. 7: (Color online) Charm quark radiative energy loss distributions S(k, x) (top panels; these are proportional to the emitted 
gluon spectra) and S(q,x) (bottom panels) for a dynamical (left) and a static medium (right), x denotes the energy fraction 
x = lo/E, and qx and kx stand for the transverse momenta |q| and |fe|, respectively. The initial quark energy is 20GeV, the 
temperature of the medium T = 225 MeV, and a path length L = 5 f m was assumed. 



such effects should improve the applicability of our approximations, by reducing large- x gluon emission. That is, the 





FIG. 8: (Color online) Same as Figure 7, but for bottom quarks. 
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total radiative energy loss should be reduced, without qualitatively affecting the energy loss ratio between static and 
dynamical media (since this ratio is seen in Fig. 6 to be largely independent of x). 

V. CONCLUSION 

Static scattering center approximation was used in all previous calculations of heavy quark radiative energy loss. 
An important consequence of this approximation is that it results in exactly zero collisional energy loss. However, 
it was recently found [11-14] that, under RHIC conditions, heavy quark collisional energy loss is significant and 
comparable to the previously calculated radiative energy loss. Since the static approximation is evidently inadequate 
in the computation of collisional energy loss, there arises a question whether such approximation is appropriate in the 
radiative energy loss case. 

We here revisited the problem of heavy quark radiative energy loss, but now in dynamical medium of thermally 
distributed massless quarks and gluons. Our work has two goals: 1) To address the applicability of static approximation 
in radiative energy loss computations, and 2) To compute collisional and radiative energy losses within a consistent 
theoretical framework. In this paper we report the first step in this direction, where we compute the 1 st order in 
opacity contribution to the radiative energy loss in a dynamical QCD medium of infinite size. 

We have shown that each individual contribution in the diagrammatic expansion of the energy loss in a dynamical 
medium is infrared divergent, due to the absence of magnetic screening [25]. However, it is interesting that the sum 
of these contributions lead to an infrared safe result. The magnetic infrared divergence is thus naturally regulated, 
eliminating the need for introducing an artificial magnetic gluon mass when computing the radiative energy loss in a 
dynamical QCD medium. 

The analytic expression for the radiative energy loss in a dynamical QCD medium was found to be remarkably 
similar to the one obtained in the static approximation. Still, the seemingly small differences, observed in the analytical 
expressions, were found to have important quantitative consequences: At the same (first) order in opacity, fast quarks 
that propagate through a dynamical QCD medium lose energy at almost twice the rate computed for a medium of 
static scatterers. Recoil of the (massless) quarks and gluons in the medium is thus a phcnomcnologically important 
effect, which can not be neglected. We found no corners of the kinematic phase-space where the static scattering 
approximation is valid, neither for light nor for heavy fast quarks. Hence, the constituents of QCD medium can not 
be approximated as static scattering centers, and for reliable predictions of radiative energy loss, dynamical effects 
have to be included. 

High precision heavy flavor measurements are expected to emerge from the upcoming high luminosity RHIC and 
LHC experiments. An important goal of heavy flavor energy loss measurements is to provide a tomographic diag- 
nostic tool for the hot QCD matter created in these collisions. Therefore, reliable quantitative predictions for these 
experiments are essential. The results presented in this paper lead to the important qualitative conclusion that the 
observed quark energy loss could be significantly larger than previously thought. Turning this qualitative insight 
into a quantitative comparison with existing data, and reliable predictions for upcoming data requires, however, a 
significant additional work. Most importantly, the present study docs not take into account coherent interference 
(LPM) effects and their modification by the finite size of the medium created in heavy- ion collision fireballs. The 
computation of heavy quark energy loss in a finite size dynamical QCD medium is therefore our next important goal. 
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APPENDIX A: ASSUMPTIONS AND APPROXIMATIONS 

1. Kinematics 

In this paper we consider a heavy quark of mass M which is produced in the remote past on its mass shell, with 
large spatial momentum p' 3> M . We choose coordinates such that the momentum of the initial quark is along the z 
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axis: 



p'=(e'*p' + ^,p>,0^ , (Al) 

We are interested in the radiative energy loss to first order in the opacity, so we study the case in which the quark 
exchanges (in arbitrary sequence) one virtual gluon with space-like momentum 

q = (go, q) = (qo, q z , q), <?o < |q| (A2) 
with a parton in the medium and radiates one (medium-modified) real gluon with time-like momentum 

k = (ko,H) = (u,k x ,k), k >\k\ (A3) 

into the medium. The quark emerges with 4- momentum p^. 

For the computation of the Feynmann diagrams given in Appendices B-D we will need cut propagators for the 
heavy quark p (D > (p)), the radiated gluon k (D> u (k)), and the exchanged gluon q (D>(q)). 

The effective 1-HTL gluon propagators for the exchanged and emitted gluons have the form given in Eq. (2.1). By 
following the procedure outlined in [25], we obtain for the cut full gluon propagator 

D> v (l) = -(l+/(/o))(iV(0Mi) + WK(0) , (A4) 

where I is gluon momentum, f(lo) = (e'°/ T — and T is the temperature of the medium. Pl.t{1) are spectral 
functions defined by 

PL .r{l) = ^ 5{f n T , L (0) - 21m ( ; 2 _n r L ( ) ^ f ) ' ( A5 ) 

It was shown in [23] that for the radiated gluon with momentum k the longitudinal contribution can be neglected 
relative to the transverse one, and that for the transverse gluon the self energy IIt(A;) can be approximated by m 2 g , 

where m g w p/^/2 = gT y/ (6+n f) / 12 is the asymptotic mass. These approximations are true in the soft rescattering 
limit u> ^> \q\ ~ \k\ ~ gT which we use in this paper. With these approximations the HTL gluon propagator for the 
emitted gluon can be simply approximated by [23] 

where is the transverse projector. The cut propagator for the radiated gluon is then given by [23, 25] 

D> v (k) « - 2n(l + f(co)) ^1 S(k - ui) , (A7) 



where w ~ \/k 2 + m 2 . 



By using Eqs. (All, A13) defined below, we obtain f(u>) — (e xE / T —l) 1 <C 1 for highly energetic jets and x > T/E. 
Eq. (A7) can then be simplified to 



D> v (k)*-2nZ^6(k -w). (A8) 
Similarly, the cut propagator for the heavy quark (with D(p) = p i_ ! ! A -i +ie ) is given by 

D>(p) = 27r^5(p Q -E). (A9) 

Unfortunately, the above approximations cannot be used for the virtual gluon mediating the collisional interaction. 
Both transverse and longitudinal contributions have to be kept in the gluon propagator D(q), and it can be shown 
numerically that both contributions are equally important. Furthermore, no further simplifications can be made in the 
expressions for the transverse and longitudinal self energies Hr(q) and Tlziq) (see Eq. (2.2)), beyond the restriction 
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(A2) to space-like momenta. However, due to this restriction, the 8 function in Eq. (A5) does not contribute to the 
exchanged gluon spectral function, leading to 

n tM . «, _ |, (1 + /M) 2Im (_^_ + gM.) . (m 

As in [33-39], we assume validity of the soft gluon (uj<^_E) and soft rescattering (|q| ~ |fe| <C k z ) approximations. 
Together with conservation of energy and momentum (p' = p + k + q) they yield 

k = (u, « *, + ^1, k z , k) , p = (£? « Pl + (fc+9 ^ +M2 , p„ -(k+qj) , (All) 

and 

/ M 2 \ 

p' = IE' ~p z +k z +q z +- ——-,p z +k z +q z ,0). (A12) 

V 2(p z +k z +q z ) ) 

In the next subsection we will show that it is reasonable to assume that q z has the same order of magnitude as |q|. 
Since |fe| <C k z and q z ~ \q\ ~ |fc|, we then also have q z ^k z . Thus k z +q z w A: z and p z +k z +q z w p z +fc 2 ~ p z +<7z ~ p z . 
Defining 

x = ^ . (A13) 



we can further rewrite 



and show that 



fc 2 + M 2 x 2 + m 2 

(p + fc) 2 - M w £ w M - (p' - fc) 2 (A14) 



P^P^(k)p v » P^P^W = P»P^{k)p' V « p'^P^ikW" « , (A15) 

where P^ v {k) is a transverse projector of radiated gluon, defined by Eq. (2.3). 
Finally, by using Eqs. (A2) and (A11)-(A12), we obtain 

fc 2 + M 2 x 2 + to 2 

£ - E -uj - q K q z - q — ^ w g z - g . (A16) 



2. qr z vs. |q| comparison 

Equation (A16) together with energy conservation implies q a <~ g z . Introducing the variable y = , f (with 
— 1 < ?/ < 1), w e can further express g z in terms of q 2 : 

The left panel in Fig. 9 shows the ratio q z /q± = q z j yfq 2 over the entire y range. We see that, except for the region 
y — > 1, q z and \q\ are comparable, and that for \y\ < 0.95 the ratio ^ remains below 3. Hence, for \y\ < 0.95, q z and 

\q\ have the same order of magnitude. 

We next want to test how important the region \y\ > 0.95 is for the energy loss. To do this, we start from the 
following equation: 

>i. (A18) 

2» V + Ben T , L (y))- + (lmn T . L (y)) 2 

which gives y-integrand of the transverse and longitudinal contributions to the energy loss (see Eqs. (Bll) and (B12)). 
By using Eq. (A18) we obtain F T ^ L {y) for a typical transverse momentum \q\ — 0.5 GeV of the exchanged gluon, which 
is shown in the right panel of Fig. 9. 

We see that the main contribution to the energy loss comes from the region \y\ < 0.95, especially after accounting 
that the contribution to the energy loss comes from the difference between longitudinal and transverse integrands (see 
Eq. (Bll) and the full curve in the right panel of Fig. 9). We additionally tested that the error made by computing the 
energy loss using our approximation in the region \y\ > 0.95 (where the approximation breaks down) is less than 2%. 
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FIG. 9: Left: The ratio = ^ as a function of y. Right: Transverse (dot-dashed) and longitudinal (dashed) contributions 
to the energy loss from virtual gluons with typical transverse momenta \q\ = 0.5 GeV, as functions of y. Full curve shows the 
difference between longitudinal and transverse contributions (see Eq. (Bll)). The vertical arrow indicates the {/-region above 
which q z exceeds \q\ by more than a factor 3. 



APPENDIX B: COMPUTATION OF DIAGRAMS M lj0 ,i - Mi, ,4 

In this appendix we present in some detail the calculation of the diagrams shown in Fig. 10. These diagrams present 
contributions where both ends of the exchanged gluon q are attached to the heavy quark, i.e. none is attached to the 
radiated gluon k and no 3-gluon vertex is involved. 

Here and later the diagrams are labeled as follows: In Mi^j, 1 denotes that these diagrams contribute to the energy 
loss to first order in opacity; i denotes how many ends of the virtual gluon q are attached to the radiated gluon k; 
and j labels the specific diagram in that class. 




I 



FIG. 10: Feynman diagrams Mi,o,i, Mi^o^, M;lo,3 and Mi,o,4 contributing to the radiative energy loss to first order in opacity. 
The large dashed circles ("blob") represent effective HTL gluon propagators [23]. A cut gluon propagator with momentum k 
and color c corresponds to the radiated gluon (ui > |k|). A cut gluon propagator with momentum q and color a corresponds to 
a collisional interaction with a parton in the medium (go < |q|). 
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1. We will first calculate the cut diagram Mf* 1 = 2ImMi j0 ,i [25]: 
M lo,i = / H g (2p'- q T)D>M(*9(2p'-qV) {p+k)2 % jp - ^ (p+fc)2 ^ _ le Hg(2fH-fc)')J?> (fc)(»g(2pf fc) g ) 
x D>(p) t a t c t cta (^f(V-M-?) ' ( B1 ) 

where D> a (k), D > (p) and D> v (q) (given by Eqs. (A8), (A9) and (A10)) arc the cut propagators for the radiated 
gluon, the heavy quark, and the exchanged gluon, respectively. Mf'g 1 then becomes 

x »d - §) (i + /(*)) (V^2to (^L. + -SgGjL) (V-.)' 

x ( .MA2 I «^2rf(p;- Po -fc„-, o )^ I g Ii . (B2) 

where we have used three of the J-functions to do the integral over (Pp. Correspondingly it should be kept in mind 
that in (B2) spatial components of p should be replaced by the corresponding components of p'—k—q which can 
be then further simplified with the approximations discussed in Appendix Al. The same will be understood when 
evaluating the other three diagrams further below. 

From Eq. (A15) and the fact that k p P pa (k) = 0we obtain 

(2p+k)PP pa {k){2p+kY « -4 ^ . (B3) 

For highly energetic jets, and by using Eqs. (2.3) and (A16), we obtain (2p'— qY P tll/ (q){2p' — q) v w 
- {2p'-qYQ^{q){2p'-qY w - 4£'V/q 2 , which leads to 

^'-" Y 2 Im (Ms + ™) * 4E " i 2 ,m (to - ™) (B4) 

By also using Eqs. (B3) and (B4), and after performing integrations over p and k , Eq. (B2) reduces to 

r 1 Ah 2 n 2 a 2 / 1 1 

Mf ,, - S > J {{p+t)2 _ M , )2 -m-p | (1 + /(„„)) 2I m ( ? - g - e - ? - 5 - ;g 

"S^ 1 "'^' (B5) 



We now use Eqs. (A14), (A16), as well as to = ^Jk 2 + m 2 g w and 4E' 2 /2E w 2£, to obtain: 

/c2 3 jfe fe 2 /" <i 4 g 

n + ff ^ f 2imn L ( g ) 2imn T ( g ) \ 

l +/W0jj. 2 l v((Z 2_ Re n L ( (? ))2 + ( Im n i (g))2 ( q 2 -ReU T ( q )) 2 + (ImU T ( q )) 2 J ' 1 ° j 
where /(go) = (e 9 °/ T — 1) _1 and T is the temperature of the medium. For small go we can expand 

1 + (B7) 
With this approximation Eq. (B6) becomes 

/d 3 k 4 k 2 

(2^ WTm^T^) 2 '< (B8) 
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where I q is given by 



A <? 2 A t\ ( 2imn L ( g ) 2imn T ( g ) 



q 7(271)4 * (?o 9z) q2 ^ + go J V(? 2 -Rcn i ( g ))2 + (Imn L ( g ))2 (g2_R e n T (g))2 + (Imn T (g))2 / 

fdqzfq_tf_(l T\( 2ImII £ (g) 2ImII T (g) \ 

~ 7 (2tt)3 g 2 +g 2 { 2 + qJ\( q 2 -ReU L ( q )) 2 + (ImU L ( q )) 2 (q 2 -ReTl T (q)) 2 + (Imn T («)) 2 J ' 1 j 

With the help of Eq. (A17) and 

= l9l77^T' ( B1 °) 

(i-y 2 ) 2 

and noting that the polarization functions Ht,l depend only on y and that 2Imn TiZ ,(y)/[(q 2 +Ren TiL (y)) 2 + 
(ImIlT,L(y)) 2 ] is an odd function of this variable, we can rewrite Eq. (B9) as 

f d 2 q 1 f 1 dy_ f 2ImU L (y) 2ImII r (y) 

q J (27r)2 27T y l V ( 9 2+Ren i ( J/ ))2 + (I m n i ( 2 ,))2 ( <7 2 +Rc n T (y))2 + (I m n T (y))2 



(Bll) 



(B12) 



rp f d*q _jl 

J (2tt)2 q2 ( q2 +/i2) ■ 

Here we used the sum rules [44] 

f 1 dyj_ 2Imn r ,L(y) 

J_ lV 2ir ( < 72 + Ren T>L ( 2 /)) 2 + (Imn TiL (y))2 \q? + R C n T , L (y=^) q 2 + Re n T , L (j/=0) 
with 

„2 

Ren TiL (y=oo) = ii- , ReII T (y=0) = 0, Re U L (y=0) = f • (B13) 
Finally Eq. (B8) becomes (2ImMi ;0 ,i = M> 01 ) 

f d 3 k d 2 q k 2 fi 2 

2ImM li0 ,i = 8Eg 4 Tt a t c t c t a J — ^—^—^ g2(g2+ ^ 2) . (B14) 

2. Next we consider the diagram M> 2 = 2ImM 1 c2 : 

M ^ = j H g (v-*n {p ,_ k)2 l _ M2+ie {pl _ k? r M2 _ ie 

xH 5 (2p+ g nD> (g)(^(2p+^ . ( B 15) 

By applying the same techniques as above and using Eq. (A14) we obtain 

/d 3 k d 2 a k 2 u 2 

(B16) 

3. Let us now compute the diagram M> 3 = 2ImMi i o,3: 

= / (-ig(2p'-qr) D> v (q) (ig(2p+ q r) (p+fc)2 _* — 2 - - (-i fl (2p+fc)") D> (fc) {igW-kY) 
- t a t c t a t c D>(p) (2n) i S^(p'-p-k-q) f fP f f« ^ 



(p'-k) 2 — M 2 — ie a c a c uvv-v - ^ ^ y; (2^) 4 (2^) 4 (2tt) 4 



- ^WcWc / (p+fc )2 j M 2 + , e {pl _ k? Z M 2 _ , £ (2 P+ ky Pp „(k)(2 P '-kn-2,) . 



2w 



x 9(1 -^)(1+ /(</„)) (2p' 

„ (5(pn — „ ,-/ / , n dp d 4 q d 4 k 

x ^^T^^-^-^^Wf W ' (B17) 
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where we have used Eqs. (A8), (A9) and (AlO). We also used three of the <5-functions to do the integral over d 3 p. By 
using Eqs. (A14)-(A16), the cut amplitude of diagram Mi j0 ,3 becomes 



/d k 4fe f d a a 

2ImII T (g) 2Imn L (q) 



,( < Z 2 -Rcn T ( <Z ))2 + (Imn T ( 9 ))2 (q 2 -RcU L (q)) 2 + (Imn t ( ? ))^ 
= -8EgH a t c t a t c Jj^_ {k2+M * x2+ml)2 I q , (B18) 

where I q is given by Eq. (Bll), giving finally 

/d 3 k d 2 a k 2 u 2 

(^R (W^ dfc)' (B19) 

4. In the same way we obtain 

2ImM li0 , 4 = -8Eg*Tt c t a t c t a J ^L- ^ __|___ q2{ ^ 2) . (B20) 

5. The sum of all four diagrams (B14), (B16), (B19), and (B20) thus becomes 

2ImMi i0 = 2ImMi i0 ,i + 2ImMi i0 ,2 + 2ImMi i0 ,3 + 2ImMi i0 ,4 

= 8Eg 4 T [t a ,t c ] [t c , t a ] J ______ q2{q2+ ^ ) , (B21) 

where [t a ,t c ] is a color commutator. 

APPENDIX C: COMPUTATION OF DIAGRAMS Mi, M - Mi,i, 4 

In this Appendix we calculate the diagrams shown in Fig. 11 where one of the ends of the exchanged gluon q is 
attached to the radiated gluon k. 

We start off with M> s l = 2ImMi i i,i: 

M> 1A = J \-ig(2p>-k'ytt) D w {k') gf cba (g pT {k'+q) x + g Xr (k- q y - g x ?{k' +ky)D> u (ig(2p+kyt c ) 

x D> a (ig(2p'- q yt a ) , , - - 2 ~* . D>( P )(2nyS^( P '- P -k-q) ■ d '" A ,H ' 



(p+k) 2 -M 2 -ie wv ' ^ y H > (2tt) 4 (2tt) 4 (2tt) 4 



where we used Eq. (A9), and performed the integral over d 4 p. As in the previous section, p = p' — k q should be 
substituted and we define G as 

G = (2p'-k'y (2p+ky (2p'- q y D w {k') D> x {k) D> T {q) (g pT {k'+ q y + g XT (k- q y - g xp (k'+ky) 

= G 1 +G 2 -G 3 (C2) 

with 

G 1 = [(2y-fc')^ w (fc')^(9)(2y-g)1[(fc'+g) A ^(fc)(2p+fc)1, 

g 2 = [(2p'-k'yD w {k') (k- q y] [(2 P +kyD> x (k)D> x °( q ) (2 P '- q y] , 

g 3 = [(2 P '-k'yD^(k')D>„ (k) (2 P +ky] [{k+k'yD>M (2 P >- q y] . (C3) 

Recalling that for the radiated gluon only the transverse part of the HTL propagator contributes, it is straightforward 
to show that (under the assumptions listed in Appendix A) the dominant contribution to Eq. (C2) comes from G3 



18 




(i.e. G\ and G2 present small corrections which can be neglected). By using Eqs. (A8), (A10) and Eq. (A16), we 
obtain 



G 3 « [^D^{k')D>Mp V ][{k+k') X D>M^P-<lY 
4p z p' z k z k' z ik ■ (k+q) 2ir 



(k' 2 +k' 2 )(k 2 + k 2 ) (k+q) 2 -m 2 +ie 2uj 



S(k - uj) 



(l + /(?o))(4fc^) 

2tt 

-IGE 2 =p- S(k - w) 



2Im(^— ^ 



"q 2 +q 2 V9 2 -n L (g) 9 2 -n T ( 9 ) 

i k-(k+q) 



2lo ~ y (k + q) 2 — to 2 + ie x 

Hence, Eq. (CI) becomes 



M> M w 8Eg 4 (if cba t b t c t 



f d*q d 
ta) J (2tt)< (27T 



(1 + /(*>)) 



1 



<7;+q 2 



21m 



1 



1 



9 2 -n L ( g ) g 2 -n T ( q ) 



(C4) 



1 



fe • (fc+q) 



x 2wS(p' -E-uj-q ) (!+/(«,)) 21m 



) 3 2w (p+fc) 2 - M 2 (fc+q) 2 -to 2 ir 

_^ „ / 1 1 

q 2 +q z 



q 2 -U L (q) q 2 -U T (q) 



Noting that if cha tbt c t a — h[t a ,t c }[t c ,t a ], the cut amplitude Mi^i then reads 



M,> M . *z?*.MK*Aj tUr* 



k 2 +M 2 X 2 + TO 2 / 



1k z 



(l+/(9o)) 



1 



fc • (fc+q) 



2ImII L (g) 



gf+q 2 
2ImIT T (q) 



(C5) 



(C6) 



' x ((p+k) 2 - M 2 ) ((k+q) 2 - TO 2 ) \(q 2 -RcU L (q)) 2 + (lmU L (q)) 2 (q 2 -RcU T (q)) 2 + (lmU T (q)Y 
The 5-function implies that 



(k+q) 2 - to 2 « -((fc+g)* + MV + to 2 ). 



(C7) 
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With the help of this and Eqs. (A14), (A16) we further obtain 
M> M « -4Eg*[t a ,t c ][t c ,t a ]J ^ ' %<>-^)(l+/U/o)) 



d 3 fc d 3 qdq _ fc ■ (fc+g) 

3 2w (2tt) 3 190 qz>[ J(q °" (k 2 +M 2 x 2 +ml)((k+q) 2 + M 2 x 2 +ml) 



Q 2 ( 2Imn L (g) 2ImII T ( g ) \ 

ql+Q 2 W-KeTl L {q)Y + (lmU L (q)) 2 (q 2 -ReU T (q)y + QmIl T {q))* J ' 

Finally, by applying the same procedure as in Eqs. (B7)-(B13), we obtain 

/d 3 k d 2 q k ■ (k+q) /i 2 

(2tt)32u> (2^)2 (fc 2 + M¥ + m 2 ) ((fc+q) 2 + M 2 x 2 + to 2 ) q 2 (q 2 +^) ' (C9) 

It is straightforward to show that the cut amplitudes of diagrams M^i^, and Mij^ each lead to the same 

result. The sum of all four diagrams computed in this section thus gives 

2ImMi ; i = 2ImMi i i ) i + 2ImMi ; i i2 + 2ImMi i i ;3 + 2ImMi ; i i4 

4 y rf 3 fc rf 2 g -2fc- (fc+g) M 2 

- 8Eg 1 [t a ,t c \[tc,ta\ J (27r)32w ( 27r )2 (fe2 + M 2 X 2 + m 2) ((fe+ q )2 + M 2 x 2 + m 2 ) q 2 (<7 2 +/x 2 ) ( j 

APPENDIX D: COMPUTATION OF DIAGRAM M 12 

In this Appendix we calculate the diagram shown in Fig. 12 where both ends of the exchanged gluon q are attached 
to the radiated gluon k: 




p' H -J"" v p< 

FIG. 12: Feynman diagram Mi,2 contributing to the radiative energy loss to first order in opacity, labeled in the same way as 
Fig. 10. 



M> 2 = J (-i 5 (2p'-A/)%)D w (fc0s/ bac ^ 

x g.f dac (g af \k'+qY + 9 af \k- q y - g a °(k'+kY)Dl u {k'){ig{2p'-qyt d ) D>(p) 
x(2nY5^( P >- P -k-q). dip ^ 



(2tt) 4 (2tt) 4 (2tt) 4 

= f bac hf dac t d / 2nS(p' -E-k -q ) H , (Dl) 

where we used Eq. (A9), and performed the integral over d 4 p. Again, p = p' k — q should be substituted and we 
define H as 

H = (2p'-fcy (2j/-fc7 D^(k')D> a (k)D> (q)DUk') 

x [g pT {k'+q) x + g XT {k- q y - g xp (k'+k) T ^j (g a0 {k' +q) a + g aP (k-q) a - g aa {k'+k)^ . (D2) 

Note that the left and right parts of the M^ 2 diagram are complex conjugates (i.e. mirror images) of each other. 
Therefore, for the three-gluon vertices on the left and right side, we go in counter-clockwise and clockwise direction, 
respectively. 
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As in the previous sections, for the radiated gluon we only consider transverse polarization. Under the assumptions 
described in Appendix A, it is straightforward to show that the dominant contribution to Eq. (D2) is given by 

h « [^D^(k')D>^k)(D^%k')r P i][(k'+krD>^ q )(k'+kr] 



.12 t/2 



V/fc' 

k' 2 ((k+q) 2 — to 2 + ie)((k+q) 2 — to 2 — it) 



2tt 



6 (ho - uj) 



2uj 



4k? (1 + /(g )) 



q 2 +q 2 



21m 



1 



1 



WE 2 



(k+q) 



g 2 -n L (<z) q 

„2 



((k+q) 2 - to 2 ) 2 



(1 + /(go)) 



q 2 +q 2 



21m 



J -n T ( q ) 
1 



q 2 -U L (q) q 2 -n T (q) 



where we have used Eqs. (A8), (A10) and Eq. (A16). 
Therefore, Eq. (Dl) becomes 



M> 2 « 8E 5 4 f bac t b f dac t d 



d 4 q d 3 k 



(2tt) 4 (2tt) 3 2w 
(1 + /(*>)) 



2tt(5(p - j B-w-9o) 



g 2 +<? 2 



21m 



((fc+g) 2 - m 2 ) 2 

By using i/ ba % = [t ,i c ], Eq. (A16) and Eq. (C7), we obtain 

d 3 k (k+q) 2 



q 2 -U L (q) q 2 -U T (q) 



/d 4 a d 
J2+Y (> 



) 3 2co ((k+q) 2 + M 2 x 2 + to?,) 2 q 2 +q : 



(!+/(*>)) 



x 27r<5(g - <7z) 



2Imn L (g) 



2Imn T (q) 



^( g 2 -Rcn L (g)) 2 + (lmU L (q)) 2 (q 2 -RcU T (q)) 2 + (Imn T (g)) 2 , 
Finally, by applying the same procedure as in Eqs. (B7)-(B13), we obtain 



r d 

2lmM 1 . 2 ^8Eg 4 T[t a ,t c }[t c t a ] / — 

J (2tt 



d 3 k d 2 q 



(k+q) 2 



) 3 2w (2tt) 2 ((k+q) 2 + M 2 x 2 + m 2 ) 2 q 2 (q 2 +^ 2 ) ' 



(D3) 



(D4) 



(D5) 



(D6) 
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